Interfacing cultured neurons to microtransducers arrays: A review of the neuro-electronic junction models

Microtransducer arrays, both metal microelectrodes and silicon-based devices, are widely used as neural interfaces to measure, extracellularly, the electrophysiological activity of excitable cells. Starting from the pioneering works at the beginning of the 70's, improvements in manufacture methods, materials, and geometrical shape have been made. Nowadays, these devices are routinely used in different experimental conditions (both in vivo and in vitro), and for several applications ranging from basic research in neuroscience to more biomedical oriented applications. However, the use of these micro-devices deeply depends on the nature of the interface (coupling) between the cell membrane and the sensitive active surface of the microtransducer. Thus, many efforts have been oriented to improve coupling conditions. Particularly, in the latest years, two innovations related to the use of carbon nanotubes as interface material and to the development of micro-structures which can be engulfed by the cell membrane have been proposed. In this work, we review what can be simulated by using simple circuital models and what happens at the interface between the sensitive active surface of the microtransducer and the neuronal membrane of in vitro neurons. We finally focus our attention on these two novel technological solutions capable to improve the coupling between neuron and micro-nano transducer

Modeling the three-dimensional connectivity of in vitro cortical ensembles coupled to Micro-Electrode Arrays

: Nowadays, in vitro three-dimensional (3D) neuronal networks are becoming a consolidated experimental model to overcome most of the intrinsic limitations of bi-dimensional (2D) assemblies. In the 3D environment, experimental evidence revealed a wider repertoire of activity patterns, characterized by a modulation of the bursting features, than the one observed in 2D cultures. However, it is not totally clear and understood what pushes the neuronal networks towards different dynamical regimes. One possible explanation could be the underlying connectivity, which could involve a larger number of neurons in a 3D rather than a 2D space and could organize following well-defined topological schemes. Driven by experimental findings, achieved by recording 3D cortical networks organized in multi-layered structures coupled to Micro-Electrode Arrays (MEAs), in the present work we developed a large-scale computational network model made up of leaky integrate-and-fire (LIF) neurons to investigate possible structural configurations able to sustain the emerging patterns of electrophysiological activity. In particular, we investigated the role of the number of layers defining a 3D assembly and the spatial distribution of the connections within and among the layers. These configurations give rise to different patterns of activity that could be compared to the ones emerging from real in vitro 3D neuronal populations. Our results suggest that the introduction of three-dimensionality induced a global reduction in both firing and bursting rates with respect to 2D models. In addition, we found that there is a minimum number of layers necessary to obtain a change in the dynamics of the network. However, the effects produced by a 3D organization of the cells is somewhat mitigated if a scale-free connectivity is implemented in either one or all the layers of the network. Finally, the best matching of the experimental data is achieved supposing a 3D connectivity organized in structured bundles of links located in different areas of the 2D network

Approximate Entropy of Spiking Series Reveals Different Dynamical States in Cortical Assemblies

Self-organized criticality theory proved that information transmission and computational performances of neural networks are optimal in critical state. By using recordings of the spontaneous activity originated by dissociated neuronal assemblies coupled to Micro-Electrode Arrays (MEAs), we tested this hypothesis using Approximate Entropy (ApEn) as a measure of complexity and information transfer. We analysed 60 min of electrophysiological activity of three neuronal cultures exhibiting either sub-critical, critical or super-critical behaviour. The firing patterns on each electrode was studied in terms of the inter-spike interval (ISI), whose complexity was quantified using ApEn. We assessed that in critical state the local complexity (measured in terms of ApEn) is larger than in sub- and super-critical conditions (mean ± std, ApEn about 0.93 ± 0.09, 0.66 ± 0.18, 0.49 ± 0.27, for the cultures in critical, sub-critical and super-critical state, respectively—differences statistically significant). Our estimations were stable when considering epochs as short as 5 min (pairwise cross- correlation of spatial distribution of mean ApEn of 94 ± 5%). These preliminary results indicate that ApEn has the potential of being a reliable and stable index to monitor local information transmission in a neuronal network during maturation. Thus, ApEn applied on ISI time series appears to be potentially useful to reflect the overall complex behaviour of the neural network, even monitoring a single specific location

Identification of excitatory-inhibitory links and network topology in large-scale neuronal assemblies from multi-electrode recordings

Functional-effective connectivity and network topology are nowadays key issues for studying brain physiological functions and pathologies. Inferring neuronal connectivity from electrophysiological recordings presents open challenges and unsolved problems. In this work, we present a cross-correlation based method for reliably estimating not only excitatory but also inhibitory links, by analyzing multi-unit spike activity from large-scale neuronal networks. The method is validated by means of realistic simulations of large-scale neuronal populations. New results related to functional connectivity estimation and network topology identification obtained by experimental electrophysiological recordings from high-density and large-scale (i.e., 4096 electrodes) microtransducer arrays coupled to in vitro neural populations are presented. Specifically, we show that: (i) functional inhibitory connections are accurately identified in in vitro cortical networks, providing that a reasonable firing rate and recording length are achieved; (ii) small-world topology, with scale-free and rich-club features are reliably obtained, on condition that a minimum number of active recording sites are available. The method and procedure can be directly extended and applied to in vivo multi-units brain activity recordings

A new simulation environment to model spontaneous and evoked activity of large-scale neuronal networks coupled to micro-electrode arrays

Introduction The use of neuronal cultures coupled to Micro-Electrode Array (MEA) is becoming a widely used and recognized experimental model for studying basic properties of information processing in neuronal systems [1]. However, the electrophysiological activity of such large-scale neuronal networks is recorded only by tens/hundreds microelectrodes. This undersampling results in a lack of information. Thus the development of a new simulation environment able to reproduce the electrophysiological behavior typically found in these preparations offers a valid help to better understand the actual dynamics. In this work, we present the main features of our software showing the simulation results of the spontaneous and evoked activity of a high-connected network

Modularity and neuronal heterogeneity: Two properties that influence in vitro neuropharmacological experiments

IntroductionThe goal of this work is to prove the relevance of the experimental model (in vitro neuronal networks in this study) when drug-delivery testing is performed.MethodsWe used dissociated cortical and hippocampal neurons coupled to Micro-Electrode Arrays (MEAs) arranged in different configurations characterized by modularity (i.e., the presence of interconnected sub-networks) and heterogeneity (i.e., the co-existence of neurons coming from brain districts). We delivered increasing concentrations of bicuculline (BIC), a neuromodulator acting on the GABAergic system, and we extracted the IC50 values (i.e., the effective concentration yielding a reduction in the response by 50%) of the mean firing rate for each configuration.ResultsWe found significant lower values of the IC50 computed for modular cortical-hippocampal ensembles than isolated cortical or hippocampal ones.DiscussionAlthough tested with a specific neuromodulator, this work aims at proving the relevance of ad hoc experimental models to perform neuropharmacological experiments to avoid errors of overestimation/underestimation leading to biased information in the characterization of the effects of a drug on neuronal networks

A single Markov-type kinetic model accounting for the macroscopic currents of all human voltage-gated sodium channel isoforms

Modelling ionic channels represents a fundamental step towards developing biologically detailed neuron models. Until recently, the voltage-gated ion channels have been mainly modelled according to the formalism introduced by the seminal works of Hodgkin and Huxley (HH). However, following the continuing achievements in the biophysical and molecular comprehension of these pore-forming transmembrane proteins, the HH formalism turned out to carry limitations and inconsistencies in reproducing the ion-channels electrophysiological behaviour. At the same time, Markov-type kinetic models have been increasingly proven to successfully replicate both the electrophysiological and biophysical features of different ion channels. However, in order to model even the finest non-conducting molecular conformational change, they are often equipped with a considerable number of states and related transitions, which make them computationally heavy and less suitable for implementation in conductance-based neurons and large networks of those. In this purely modelling study we develop a Markov-type kinetic model for all human voltage-gated sodium channels (VGSCs). The model framework is detailed, unifying (i.e., it accounts for all ion-channel isoforms) and computationally efficient (i.e. with a minimal set of states and transitions). The electrophysiological data to be modelled are gathered from previously published studies on whole-cell patch-clamp experiments in mammalian cell lines heterologously expressing the human VGSC subtypes (from NaV1.1 to NaV1.9). By adopting a minimum sequence of states, and using the same state diagram for all the distinct isoforms, the model ensures the lightest computational load when used in neuron models and neural networks of increasing complexity. The transitions between the states are described by original ordinary differential equations, which represent the rate of the state transitions as a function of voltage (i.e., membrane potential). The kinetic model, developed in the NEURON simulation environment, appears to be the simplest and most parsimonious way for a detailed phenomenological description of the human VGSCs electrophysiological behaviour

Stimulus-Evoked Activity Modulation of In Vitro Engineered Cortical and Hippocampal Networks

The delivery of electrical stimuli is crucial to shape the electrophysiological activity of neuronal populations and to appreciate the response of the different brain circuits involved. In the present work, we used dissociated cortical and hippocampal networks coupled to Micro-Electrode Arrays (MEAs) to investigate the features of their evoked response when a low-frequency (0.2 Hz) electrical stimulation protocol is delivered. In particular, cortical and hippocampal neurons were topologically organized to recreate interconnected sub-populations with a polydimethylsiloxane (PDMS) mask, which guaranteed the segregation of the cell bodies and the connections among the sub-regions through microchannels. We found that cortical assemblies were more reactive than hippocampal ones. Despite both configurations exhibiting a fast (<35 ms) response, this did not uniformly distribute over the MEA in the hippocampal networks. Moreover, the propagation of the stimuli-evoked activity within the networks showed a late (35-500 ms) response only in the cortical assemblies. The achieved results suggest the importance of the neuronal target when electrical stimulation experiments are performed. Not all neuronal types display the same response, and in light of transferring stimulation protocols to in vivo applications, it becomes fundamental to design realistic in vitro brain-on-a-chip devices to investigate the dynamical properties of complex neuronal circuits

In Vitro Studies of Neuronal Networks and Synaptic Plasticity in Invertebrates and in Mammals Using Multielectrode Arrays

Brain functions are strictly dependent on neural connections formed during development and modified during life. The cellular and molecular mechanisms underlying synaptogenesis and plastic changes involved in learning and memory have been analyzed in detail in simple animals such as invertebrates and in circuits of mammalian brains mainly by intracellular recordings of neuronal activity. In the last decades, the evolution of techniques such as microelectrode arrays (MEAs) that allow simultaneous, long-lasting, noninvasive, extracellular recordings from a large number of neurons has proven very useful to study long-term processes in neuronal networks in vivo and in vitro. In this work, we start off by briefly reviewing the microelectrode array technology and the optimization of the coupling between neurons and microtransducers to detect subthreshold synaptic signals. Then, we report MEA studies of circuit formation and activity in invertebrate models such as Lymnaea, Aplysia, and Helix. In the following sections, we analyze plasticity and connectivity in cultures of mammalian dissociated neurons, focusing on spontaneous activity and electrical stimulation. We conclude by discussing plasticity in closed-loop experiments

An FPGA platform for real-time simulation of spiking neuronal networks

In the last years, the idea to dynamically interface biological neurons with artificial ones has become more and more urgent. The reason is essentially due to the design of innovative neuroprostheses where biological cell assemblies of the brain can be substituted by artificial ones. For closed-loop experiments with biological neuronal networks interfaced with in silico modeled networks, several technological challenges need to be faced, from the low-level interfacing between the living tissue and the computational model to the implementation of the latter in a suitable form for real-time processing. Field programmable gate arrays (FPGAs) can improve flexibility when simple neuronal models are required, obtaining good accuracy, real-time performance, and the possibility to create a hybrid system without any custom hardware, just programming the hardware to achieve the required functionality. In this paper, this possibility is explored presenting a modular and efficient FPGA design of an in silico spiking neural network exploiting the Izhikevich model. The proposed system, prototypically implemented on a Xilinx Virtex 6 device, is able to simulate a fully connected network counting up to 1,440 neurons, in real-time, at a sampling rate of 10 kHz, which is reasonable for small to medium scale extra-cellular closed-loop experiments